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We investigate the possibilities of preserving and manipu- 
lating the coherence of atomic two level systems by ideal pro- 
jective homodyne detection and feedback. For this purpose, 
the photon emission process is described on timescales much 
shorter than the lifetime of the excited state using a model 
based on Wigner-Weisskopf theory. The back-action of this 
emission process is analytically described as a quantum diffu- 
sion of the Bloch vector. It is shown that the evolution of the 
atomic wavefunction can be controlled completely using the 
results of homodyne detection. This allows the stabilization of 
a known quantum state or the creation of coherent states by a 
feedback mechanism. However, the feedback mechanism can 
never compensate the dissipative effects of quantum fluctua- 
tions even though the coherent state of the system is known 
at all times. 

PACS numbers: 42.50.Lc, 03.65.-w 



I. INTRODUCTION 

Homodyne detection is a procedure which allows the 
measurement of a quadrature component of the light 
field. In quantum optics, it has long been applied to in- 
vestigate squeezing effects of non-linear optical systems. 
In that context it is not necessary to understand the time 
resolved properties of the homodyne detection process it- 
self, since the measurement is a time average over many 
photon emissions from the light field source under inves- 
tigation. If homodyne detection is to be applied to the 
observation of individual quantum systems, however, a 
time resolved quantum mechanical treatment of the pro- 
jective measurements performed by a homodyne detec- 
tion setup is necessary. Such a description not only helps 
to improve our understanding of the quantum mechanics 
of photon emissions, but also provides new methods of 
controlling the dynamics of quantum systems. 

A complete quantum theory of homodyne detection 
must take into account both the quantization of the field 
modes and the temporal evolution of the continuous light 
field entering the detector. The problem of quantization 
has been addressed in a number of publications [Q-l). 
An approach to the problem of measuring a system as it 
evolves in time is presented by H. Carmichael in Q. It 
has been applied to a number of measurement scenarios, 
including homodyne and heterodyne detection as 
well as to feedback scenarios Pi- 



In section ^ of this paper, we present an alternative 
derivation of this quantum trajectory approach to time 
resolved observation of quantum systems by homodyne 
detection, starting from the projective measurement pos- 
tulate formulated in In section [II this measure- 
ment theory is then applied to a description of the spon- 
taneous emission process based on Wigner-Weisskopf the- 
ory, which takes into account the field dynamics of the 
electromagnetic continuum of modes. This derivation 
provides a shortcut to the formulation of the back-action 
of the many photon measurements on the two level atom 
and illustrates the problems of time resolution and of 
non-locality in quantum trajectories, which are discussed 
in Q . Once this control over the informatin loss induced 
by spontaneous emission has been established, the ob- 
served changes in the atomic wavefunction can be com- 
pensated by feedback. This is investigated in section IV. 
In section |y|, the results are interpreted and conclusions 
are presented. 



II. A QUANTUM MECHANICAL DESCRIPTION 
OF TIME RESOLVED HOMODYNE DETECTION 

A. Time resolution and projective measurements in 
the homodyne detection setup 

If we know that the lincwidth of the emission we wish 
to observe is much smaller than 1/t, then we can con- 
clude that a time resolution of r will be sufficient to 
observe the dynamics of photon emission from the sys- 
tem. We therefore assume that the detectors used pro- 
vide us with reliable information about the photon num- 
bers which entered each of the detectors during a time 
interval t. In fact this is much closer to the experimental 
situation than the assumption of a continuous measure- 
ment in which the arrival of a photon is detected with infi- 
nite time resolution. For all time resolved measurements 
there is a time window of length r which will usually not 
be shorter than a femtosecond. We can then describe 
the homodyne detection as a projective measurement on 
a system of two light field modes of frequency luq. The 
state to be measured is the product state of the input 
field state | 0) and the coherent state | a) produced by 
the local oscillator, a is the complex field amplitude of 
the local oscillator mode emitted during the time interval 
T. This amplitude is related to the intensity / emitted by 
the local oscillater via the relation / = a*a/r. The mea- 
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surement base is the photon number base of the symmet- 
ric and anti-symmetric hnear combinations of the input 
mode and the local oscillator mode, | n+,n_). In terms 
of the input and local oscillator fields, this state may be 
written as 



it + 6t 



St 



V2 



vacuum) 



(1) 



where a) is the creation operator of the local oscillator 



light field mode {{a\a 



t - 



(a I) and 6^ is the creation 



operator of the input field. 

The effect of this measurement on the input field state 
I (^) is a projective measurement on a set of states which 
is neither orthogonal nor normalized within the subspace 
of interest. The state | P(ri+,rt_)) projected into by a 
measurement of n+ and ri_ photons in the detectors is 
found by forming the scalar product with {a \ in the 
subspace of the local oscillator field. Using the fact that 
(a I is an Eigenvector of a\ one obtains 

I P(n+,n„)) = 



5M / a* - fet 



V2 



V2 



vacuum) 



(2) 



Note that even though this state is not normalized, the 
probability of a measurement of n+ and n_ is neverthe- 
less given by 



(3) 



The projection operators associated with these states 
form a complete probability measure for the subsystem 
of the input light field. It has been shown that the pho- 
ton number difference n+ — n_ corresponds very well to 
the quadrature component of the input field which is in 
phase with the local oscillator if the input field photon 
number is much smaller than a* a (see for example [^|). 
Therefore, the classical interpretation of a homodyne de- 
tection as a measurement of a quadrature component of 
the input field also applies to the quantum mechanical 
measurement. 



B. Application to low intensity light fields 

The light field sources we wish to investigate in the 
following will usually emit only zero or one photon (ni„ = 
0, 1) during the time interval r. If the input state is the 
vacuum state, the measurement probability is given by 

) = I {vacuum \ P(n+,n_)) | 



n-(-!n_! 



a a 



(ri++n-) 



(4) 



For a fixed sum of photons, iV = n_|_ -I- n_, this probabil- 
ity distribution is the binomial distribution which results 
from the random scattering of the photons by the beam 
splitter. For large photon numbers ^ 1, this proba- 
bility distribution may be approximated by a Gaussian. 
The probability of measuring the photon number differ- 
ence An = n_|_— n_ regardless of the total photon number 
is then given by 



Pvac 



.(An) 



1 



: exp[- 



2a*a' 



(5) 



This probability distribution represents the vacuum fluc- 
tuations of the in-phase quadrature of the incoming field. 
In the limit of very strong local oscillator fields, the mea- 
sured value of this quadrature corresponding to the pho- 
ton number difference of An is given by An/2 | a |, 
reproducing the expected quantum uncertainty of 1 /4 in 
the quadrature component. 

In the case of weak one photon contributions, 



$/3) ~\ vacuum) + (3 \ni. 



I) 



(6) 



the probability distribution is modified only slightly by 
the contributions linear in (3. To simplify the formalism, 
the non-orthogonal measurement base of the homodyne 
detection may be reduced to the components lying in the 
two dimensional subspace of zero or one photon and the 
amplitude factor corresponding to the binomial distribu- 
tion for the vacuum can be approximated by a Gaus- 
sian. The resulting measurement base depends only on 
the photon number difference An and is given by 

I P(An)) = (27ra*a)-i/4 | P(An)) 



(27ra*a)-^/*exp[ 
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4a* a 
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(7) 



Note that the weak field condition can be fulfilled for 
any input field intensity by chosing a time scale r which 
is much smaller than the average rate of photons corre- 
sponding to the intensity of the input field. Therefore, 
the simplified projection state | P(An)) can describe the 
time evolution of any homodyne detection scenario. Ap- 
plying equation (0) to the weak field state | results 
in a probability distribution of 



Pf}{An) = 



1 



^/2^ 



:exp[ 



{An - {a*f3 + P*a)y 
2a*a 



(8) 



Effectively, the weak coherent field characterized by /3 
shifts the Gaussian distribution by just the amount ex- 
pected from classical homodyne detection. 
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III. QUANTUM DIFFUSION OF A TWO LEVEL 
ATOM 

A. Description of the emission process 

The coherent quantum dynamical evolution of an 
atomic system interacting with the light field continuum 
is given by Wigner Weisskopf theory. The spatio- tem- 
poral interpretation of this theory shows that emission 
processes can be described by a temporal evolution of 
the following type : 

{E- vacuum \ = ei-i'^o-T/2)t ^g^^ 

I otherwise 

(9b) 

where | E\ vacuum) is the state of the excited atom in the 
light field vacuum and | G; r) is the state of the ground 
state atom with a photon at a distance of r from the 
atom, r is the rate of spontaneous emission. 

During a time interval r which is much shorter than 
1/r, the product state of the light field vacuum and an 
arbitrary linear combination of the excited state | E) and 
the ground state | G), 

I ^(0)) — ce \ E] vacuum) + cg\ G; vacuum), (10) 

therefore evolves into an entangled state of the system 
and the light field, with the light field being in the vac- 
uum state or in a single photon state. The single photon 
wavefunction will be located within a distance of ct from 
the atom. Since the emission has a well defined frequency 
and angular dependence, the mode into which the pho- 
ton has been emitted is a well defined mode and the total 
state of atom and field can be written as 

I *(t)) = Ce{1 - {iuJo + r/2)T) I E; vacuum) 

+cg I G; vacuum) + ceVTt | G; no = 1). (11) 

The amplitude factor of VTV is found by normalizing the 
rectangular mode emitted by the system after a time in- 
terval of T. This is done by dividing the amplitude of 
\fvjc emitted into the spatially continuous one dimen- 
sional field by the amplitude of the normalized rectangu- 
lar mode of length ct which is equal to \f\Jcr. Note that 
the amplitude of the emission depends on the square root 
of T, reflecting the linear increase of emission probability 
with time. 

In homodyne detection, the frequency of the local os- 
cillator is also loq. Therefore, it is useful to transform to 
the interaction picture using the time-dependent trans- 
formation 

\E) =c-'^''^ \E). (12) 



This effectively removes the terms oscillating with wq 
from the system dynamics by describing the phase rela- 
tion between the excited state component and the ground 
state component of the system state not in terms of an 
absolute phase but in terms of the phase relative to the 
local oscillator. Note that heterodyne detection may also 
be described by this formalism if the time scale r is cho- 
sen so that the detuning between the local oscillator 
and the system dynamics satisfies the requirement that 
8ijJT ^ 1. The dynamical evolution of the system phase 
relative to the local oscillator phase can then be included 
in the measurement scenario. 

The equation for the evolution of the wavefunction dur- 
ing short time intervals r is 

I ^(j)) = C£;(l - rT/2) I E; vacuum) 

-|-cg I G; vacuum) + ceVTt \ G; no = 1). (13) 

The projection postulate of homodyne detection can now 
be applied to the light field part of this correlated system- 
field state, resulting in an effective projection of the state 
of the atomic system into a state which depends on the 
measurement result of the homodyne detection. 

Although the real physical measurement will take place 
many time intervals t later, after the light field wavefunc- 
tion has traveled the distance to the detectors, it is pos- 
sible to interpret the measurement as an instantaneous 
projection on the local light field state, since the light 
field signal, once emitted, will not interact with the atom 
again. This artificial choice of the instant in which the 
state is projected is consistent with the basic theory of 
quantum measurement, since it is not possible to distin- 
guish between a projection at the time of a measurement 
and a projection which anticipates the measurement. It 
is only our subjective expectation of causality which leads 
us to prefer placing the collapse after the measurement. 

Also the measurement described here is not continu- 
ous in the sense discussed in |6|j^ . The discrete measure- 
ments performed are measurements of properties of the 
total interval t. For this reason r is not written as dt, 
which would suggest the limit of infinite time resolution. 
The experimental result corresponding to the situation 
described here is a series of photon numbers without any 
zero photon time intervals separating the time windows 
of each measurement. 

If the wavefunction is regarded as an epistemological 
tool describing not physical reality itself but only our 
knowledge of it, then the projective measurement at the 
atom is simply an expression of the information gained 
about the events at the atom, independent of the time 
at which the information is actually obtained. In a com- 
pletely relativistic theory of measurement, this episte- 
mological effect of future knowledge on the state of the 
past must be considered, as was already pointed out by 
Einstein ||ll|. The arbitrary subdivision of the flow of 
time into segments of duration r is also an epistemolog- 
ical consequence as it represents the time aspect of the 
space-time measurement base defined by the experimen- 
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tal setup. Information about the atomic system dynam- 
ics on shorter time scales is not obtained and cannot be 
included in the evolution of the system wavefunction. 



B. Influence of the homodyne detection on the 
system dynamics 

For an arbitrary system with known initial wavefunc- 
tion, the measurement protocol of repeated homodyne 
detections provides a complete description of the evo- 
lution of the system wavefunction. In this sense, the 
method described here is a generalized quantum trajec- 
tory approach. However, instead of using a master equa- 
tion approach, Schroedingers equation of the system-field 
interaction is solved and the homodyne detection events 
are described as projective measurements on the corre- 
lated system-field wavefunction. The detector simulated 
by this approach is not a single photon counting device 
but rather corresponds to a photo diode with a high 
time-resolution. Indeed, the measurement need not re- 
solve single photon counts to achieve a useful precision 
for predictions about the quantum system. 

To investigate the dynamics of the atomic wavefunc- 
tion observed by homodyne detection, it is useful to ex- 
amine the effect of a single projective measurement of 
photon number difference An on an arbitrary system 
wavefunction after an emission time segment of r. Before 
the projective measurement, the correlated system-field 
wavefunction is given by equation (p^. After the pro- 
jective measurement of An the system is in a pure state 
described in the interaction picture by 



^(r)) = (P(An) I ^{r)) 



(27ra*a)-^/Sxp[ 



An" 

4a* a 



An, 



X i^CEil - rr/2) I E) + (CG + ceVTt — ) \ G) j . (14) 

The squared length of this state vector is a measure of 
the probability of measuring (n^ , n_ ) . Proceding as in 
the derivation of equation (||) and assuming that Fr ^ 1 , 
it is possible to represent the deviation of this probabil- 
ity distribution from the vacuum state distribution by a 
Gaussian distribution shifted by a term linear in vTV: 

p{An) w 



1 



\/27ra*a 



cxp[ 



(An- I a I Vt^{cec*g + c^cg))" 
2a*a 



(15) 



Here and in the following, the phase relation between the 
local oscillator and the atomic dipole is defined so that it 
is zero if both ce and Cg are real and positive. The small 



deviation from the quantum vacuum case are clearly re- 
lated to the dipole expectation value of the atomic sys- 
tem. Consequently, they indicate the dipole field emitted 
by the system. The normalized change in the state of the 
system is orthogonal to the initial state | "0(0)) for small 
changes. It is found by projecting | iP{t)) on the sub- 
space orthogonal to | V'(0))j normalizing by dividing by 
the amplitude of the parallel component. 



V>(t))- |0(O))(V>(O) |0(t)) 
(V'(O) I 0(r)) 



(16) 



The normalized expression for the system wavefunction is 
then given by | ^(0))-|- | 5tp{T)). Applying the condition 
that Ft <C 1 we neglect all terms above second order in 
"v/Ft. The change of the wavefunction | S^{t)) within 
the time interval r conditioned by the measurement of 
An is then given by 



I S^{r)) = -Vf^^cl (c*a I E) 
\a\ V 



+Fr— ^c^; (c^ I E) 



ce\G) 



^CECG {c*G \E)-C*e\ G)) . (17) 

The dominant effect within one time interval is given 
by the diffusion term proportional to ^/TrAn. The av- 
erage of these diffusion steps is only proportional to Fr 
however, since the expectation value of An /a is of the or- 
der of vTV. On the measurement timescale r, the quan- 
tum fluctuations therefore cause a random walk type dif- 
fusion of the atomic state, just as one would expect from 
classical noise. However, the length and the phase of the 
diffusion step are functions of the initial system state. 
On timescales of 1/F, the non-zero average of An/a, the 
diffusion term proportional to the square of the measure- 
ment, and the deterministic drift terms contribute to the 
change in the system state. All of these terms give rise 
to a slow drift of the wavefunction towards the ground 
state. To illustrate the dependence of the diffusion con- 
stant on the state of the atomic system and to identify 
the drift terms we will now formulate this dependence of 
these processes in terms of the Bloch vector. 



C. Quantum diffusion of the Bloch vector 

The dynamics of two level systems can be visualized 
using the Bloch vector representation. This three dimen- 
sional vector incorporates the excitation (population in- 
version), the dipole and the dipole current of the two 
level atom as its orthogonal components. When the time 
dependent transformation given in equation (|l|) is used, 
the two components orthogonal to the excitation Sz de- 
scribe the in-phase and the (tt/2) out-of-phase compo- 
nents, respectively, of the dipole oscillations relative to 
the local oscillator. 
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= 2Re{{i;\E){G\^)) 
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(18a) 
(18b) 
(18c) 



If a is a real number, Sx is the in-phase component of 
the atomic dipole and Sy is the out-of-phase component. 
The diffusion step associated with a measurement resuh 
of An in terms of the Bloch vector is derived from | ■(/'(O)) 
and I 6'iP{t)) by using 

5s, = 2i?e((V(0) \E){G\ 8^{r)) 

+ {5^{T)\E){G\m)) (19a) 
5sy = 2/to((?/'(0) \E){G\ S^{t)) 

+ {S^/j{T)\E){G\m)) (19b) 
5s, = 2i?e((V'(0) \E){E\ 57A(t)) 

-{m\G){G\S,p{r))). (19c) 

The change in the Bloch vector 6s can be expressed 
in terms of the Bloch vector s corresponding to | "0(0)). 
The result is the diffusion step of the Bloch vector, 



5s^ 

SSy 

6s. 




(20) 



As in the case of the state vector representation, the 
Bloch vector diffusion step is composed of a fast random 
walk diffusion on the time scale of r and drift terms on 
a time scale of l/F. What is more apparent in the Bloch 
vector representation is the possible separation into phase 
sensitive and phase independent changes. 



D. Interpretation of the contributions to the 
diffusion step 

To interpret the diffusion step, it is useful to analyze 
the seperate contributions in more detail. In particular 
it is helpful to investigate the diffusion step component 
representing the random walk caused by the fluctuating 
light field. 



6sx 



An 



1 



fluctuation 



SxSy 

-Sx(l + Sz) 



(21) 



Figure |T| illustrates this diffusion on the Bloch sphere. In 
order to analyse the diffusive motion of the Bloch vector, 
one may separate the absolute value from the direction 
of the diffusion. The diffusion constant is then given by 



(22) 



The magnitude of the diffusion is therefore independent 
of the phase relation between the local oscillator and the 
atomic dipole. It does depend on the excitation of the 
atom however. It is maximal for the excited state and 
zero for the ground state, indicating that the ground state 
will not interact with the field vacuum. The direction of 
6s is given by 



6s 



6s 




(23a) 



where ej^ is the unit vector perpendicular to both s and 
the z-axis and e|| is the unit vector perpendicular to both 
s and e_L. 



s^ + 



and 



4 + 4 



SySz 



(23b) 



As can be seen from these results, a Bloch vector which is 
in phase with the local oscillator will show random rota- 
tions in the Sx, s, plane. A Bloch vector polarized out of 
phase with the local oscillator by ±7r/2 {sx — 0) will un- 
dergo phase diffusion only, with s, remaining constant. 
The ratio between the diffusion constant of pure phase 
diffusion and the diffusion constant of excitation diffusion 
is given by (sy/sx)'^ — tan((/))^, where (f) is the phase dif- 
ference between the local oscillator and the dipole oscilla- 
tions of the atom. Note that a randomly varying classical 
in phase field would give rise to Rabi oscillations around 
Sy with a step length proportional to -\/l — s^. The prop- 
erties that only phase diffusion occurs for Sx — and that 
the rotations in the Sy — plane preserve phase are also 
properties of Rabi rotations induced by classical fields. 
However, the step length dependence of 1 -I- s, clearly in- 
dicates a difference between the effects of quantum noise 
and of classical noise. 

On a timescale of l/F, the contribution of the random 
walk is given by the non-vanishing expectation value of 
the homodyne detection photon number difference. The 
probability distribution of the measurement results is ap- 
proximately given by a Gaussian with 
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( I — r> = VrTsa; 
a 



= 1. 



(24a) 



(24b) 



On a timescale of l/F a large number of measurements 
will have been performed since r <C 1/r. Therefore, the 
net change of the system state can be evaluated using 
these expectation values. The random walk drift is then 
given by 



{SSy) 



fluctuation 



TtSx I -s.^Sy I . (25) 



This is exactly compensated by the other phase sensitive 
contribution effective on a time scale of l/F, 



{SSy) 
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(26) 



This compensation has a clear physical reason. The dy- 
namics of the system should only depend on the incoming 
field. However, the measurement is being performed on 
the field coming from the atomic system which is a sum 
of the fields passing the atom and the dipole field emit- 
ted by the atom. Although it is impossible to seperate 
the contributions in a quantum mechanical measurement, 
the long term averages can actually be corrected by sub- 
tracting the averaged dipole field contribution from the 
measured fields. In effect, this leaves only two drift terms 
which infiuence the l/F time scale dynamics. 

The drift term which is not dependent on the measure- 
ment result Art is given by 



„2 



(27) 



dipole 



This term describes a reduction of Sz as a function of 
si + Sy. Since Sz is the energy expectation value of the 
atomic system and -I- Sy is the square of the atomic 
dipole, this process describes a loss of energy correspond- 
ing to the emission of radiation from a classical oscillating 



dipole. It is therefore eq 
ing super radiance [|12| 



: uivalent to the equations describ- 

jjfsj, as well as to a semiclassical 

textbook approach to spontaneous emission jljl . Unlike 
the other terms, this term is entirely free of quantum 
fluctuation effects. It is completely deterministic and de- 
pends only on the expectation values of the atomic sys- 
tem state. It is a fascinating feature of the quantum 
theoretical formulation of spontaneous emission applied 
in this investigation that it automatically produces such 
a semiclassical term from the simple photon emission de- 
scribed by equation (0). 



The three terms discussed up to this point can be un- 
derstood in terms of the action of the incoming quantum 
fluctuations on the system and in terms of the semiclassi- 
cal emission caused by the dipole oscillations. The fourth 
term is more difficult to understand since it depends on 
the quantum fluctuation measurements but is not sensi- 
tive to the phase of the local oscillator. It seems to be 
a mixed effect of dipole emission and quantum fluctu- 
ations, possibly related to the quantum fluctuations of 
the atomic dipole. The contribution of this mixed term 
to the dynamics is given by 




Fr , An2 



>(1 



SzSy 
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mixed 



(28) 



The sum of the drift dynamics to be expected on a 
timescale of l/F is given by the sum of the dipole con- 
tribution and the mixed term. The total drift vectors 
resulting from this sum arc shown in figure 0. 



E. Contributions to the exponential decay of the 
average excitation 

At first it may seem confusing that a simple sponta- 
neous emission from an isolated atom should give rise to 
such a complicated variety of dynamical effects. After 
all, the Bloch vector dynamics of an ensemble of atoms 
is described by 




(29) 



If one separates the time derivative of the Bloch vector 
into a component orthogonal to the Bloch vector and one 
parallel to it, however, the connection to the dynamics 
observed in the homodyne detection becomes apparent: 





(30) 



For I (s) 1=1, the expected reduction of the length of the 
Bloch vector is an effect of the diffusion caused by the 
Ss fluctuations term. Its rate is given by one half of the 
diffusion constant F (1 -t- Sz)'^- The rotation of the Bloch 
vector is given by the sum of the two drift terms, Ssdipoie 
and Ssmij;f,(i- For | (s) |= 1, the average expected change 
in Sz can be separated into contributions related to the 
diffusion of the Bloch vector, to the dipole emission, and 
to the mixed term, 
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(31a) 

(31b) 
(31c) 
(31d) 



Figure ^ shows this change in the relative contributions 
to the exponential decay of the atomic excitation as a 
function of Sz- For the excited state, Sz — +1, the ex- 
ponential relaxation of Sz is a result of the diffusion of 
the Bloch vector due to quantum fluctuations, while for 
states close to the ground state, Sz ~ — 1 , the exponential 
relaxation to = — 1 is dominated by the semiclassical 
dipole emission. For maximally polarized states, Sz = 0, 
the emission contributions of the dipole term and of the 
mixed term are equally strong, while the diffusion has no 
effect. Note that for — 1 < < 0, the diffusion caused 
by the quantum fluctuations actually tends to excite the 
atom. However, this effect is compensated by the mixed 
terms and the dipole emission. 



IV. COMPENSATING QUANTUM 
FLUCTUATIONS BY FEEDBACK 



A. Feedback setup 



Section [II describes the measurement of quantum fluc- 
tuations propagating away from an atomic two level sys- 
tem by homodyne detection and the implications of the 
measurement for the dynamics of the atomic state. Since 
maximal knowledge is obtained about the field state, the 
atomic system state is in a pure state after the mea- 
surement. This suggests the possibility of using feed- 
back based on the measurement results as a means to 
manipulate the quantum state. A general formalism for 
this procedure based on the master equation approach to 
quantum trajectories has been presented in |p| and was 
applied to light fields in an open cavity in H. In the 
following, we will discuss the possibility and the physical 
problems involved in applying a feedback scheme to the 
two level atom. 

On a timescale of r, the homodyne detection does not 
give any information about the state of the atom. p{An) 
is nearly equal to the vacuum probability distribution. 
The information obtained in a single field measurement 
is therefore not information about the atomic system it- 
self. Instead, it is information about the quantum fluc- 
tuations acting on the system. This sensitivity to the 
vacuum fluctuations (as opposed to the system state) is 
ideal for stabilizing quantum states by reversing the ef- 
fects of the measurement process [09] . If we were dealing 



with a classical system it would be possible to measure 
the exact forces involved. A feedback mechanism could 
then compensate the forces, thereby removing the effects 
of the incoming noise. In quantum mechanics however, 
the quadratures of the light field do not commute. There- 
fore the force acting on the atomic system can never be 
known completely. In the balanced homodyne detection 
setup discussed here we only know the quadrature of the 
light field which is in phase with the local oscillator. In 
a classical system, this lack of information would make 
control impossible. In the quantum system however, we 
can still achieve perfect control of the atomic wavefunc- 
tion. 

In order to manipulate the known atomic state the 
electromagnetic field at the atom may be modified. For 
example, a resonant coherent driving field may be cou- 
pled to the atomic system. In order to ensure that this 
field has a fixed phase relation with the local oscillator 
used in the measurement it would be natural to utilize 
the same coherent light field source in the homodyne de- 
tection setup and in the feedback loop. Since the local 
oscillator field is very strong, only a negligibly small por- 
tion needs to be redirected in order to achieve control of 
the system dynamics. In fact, since the strength of the 
field expectation value needs to compensate for quantum 
fiuctuations only, the portion of the local oscillator in- 
tensity used for this purpose is approximately given by 
l/(4a*Q!). Note that it is possible to add the effect of 
the Rabi rotations caused by the feedback field and the 
effects of the quantum fluctuations because the Heisen- 
berg equations of motion for the fully quantum mechani- 
cal field-atom interaction are linear in the field variables. 
The dynamical effects of the feedback and of the diffusive 
evolution are therefore seperable. 

A controllable reflector can be used to coherently ma- 
nipulate the system depending on the measurement re- 
sults of the homodyne detection. This feedback modifies 
the dissipative dynamics of the atomic system. If the de- 
lay between emission and feedback is much shorter than 
1/r, it is possible to compensate the effects of fluctua- 
tions on a known system state. Note that this either re- 
quires a low decay rate F or a very fast feedback loop. In 
an optical setup a typical unmodified lifetime of nanosec- 
onds would require a feedback loop much shorter than 10 
cm in length, so the light field signal can return in time, 
with the purely optical dissipative photon detectors inte- 
grated into this loop. 

The effect of a delay time of At — rriT between the 
emission and the arrival of the corresponding feedback 
signal at the atom can be estimated by considering the 
number of uncompensated diffusion steps m which oc- 
cur during the delay time. The probability distribution 
over the sum of the measurement results Sn accumulated 
within that uncontrolled interval is given by a Gaussian 
with a standard deviation of | Qf | . This may be mul- 
tiplied with the expected uncompensated diffusion step 
length at the stabilized point on the Bloch sphere in or- 
der to determine the reliability of feedback stabilization 
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with a non vanishing delay time. In the following, how- 
ever, we shall concentrate on the description of effectively 
instantaneous feedback. 

Because the total light field propagating away from the 
atomic system is measured in the homodyne detection, 
the reflected control field will also be measured in the 
subsequent time interval. This must be subtracted from 
the measurement signal for further feedback, since the 
feedback effect itself should not be compensated. Note 
that the linearity of the field dynamics involved actually 
allows this seperation despite the fact that the feedback 
field and the quantum fluctuations interact with the sys- 
tem in a qualitatively different manner as described be- 
low. Intuitively, one would expect a complete compensa- 
tion of the fluctuations if the averaged measurements of 
An are zero. However, we will see that this is not so, ow- 
ing to the fact that only one quadrature of the quantum 
fluctuations can be compensated. 



For Sy = 0, this condition is fuUfiUed if 



/(An) = -(l + s,) 



An 

2 I a 



(35) 



where is the Bloch vector component of the stabilized 
state. Note that /(An) is not a function of the present 
state of the atomic system. 

A particularily interesting case is obtained when the 
feedback is chosen to compensate for the maximal fluctu- 
ation effects possible, stabilizing the excited state of the 
two level atom. If the system is not in the excited state, 
the diffusion step associated with a homodyne detection 
event of An is now given by the sum of the original dif- 
fusion step and the feedback. 




At 



a 



eff. 



B. Effects of feedback on timescales shorter than l/F 

Although a number of feedback scenarios may be dis- 
cussed even in the simple context of the two level atom, 
we will now concentrate on the possibility of stabilizing 
a known state with Sj, = by effectively instantaneous 
feedback. The feedback is given by a coherent driving 
field inducing Rabi rotations around the Sy axis. The 
effect of this feedback experienced by the system during 
the time interval r can be written as 



2VrT/(An) 



feedback 







(32) 



where /(An) is the feedback field describing the response 
of the coherent control to the most recent measurement 
result. An is the measurement result obtained from the 
quantum fluctuations acting on the system in the previ- 
ous time interval. Since the time intervals are small on 
the scale of the system dynamics, however, we will sim- 
plify matters by summing the effects of the diffusion step 
and its subsequent feedback as if they occurred in the 
same time interval to obtain the effective total diffusion. 

The homodyne detection measurement in the following 
time interval will be modified by the feedback field. The 
change in the An measurement caused by the feedback in 
the next measurement, Snext, uiay be determined using 
equation(||): 



/(An). 



(33) 



In order to stabilize a given system state s, /(An) must 
be chosen to compensate the diffusion step, i.e. 



fluctuation 




0. 



(34) 




An 



SxSy 



(36) 



As shown in figure ^, this effective diffusion step is the 
exact inversion of the diffusion step without feedback. 
The roles of Sz = —1 and of Sz = +1 have been re- 
versed. Consequently, the sensitivity of the system to 
quantum fluctuations is now maximal at the ground 
state, Sz = —1. The diffusion constant is proportional 
to (1 — Sz)'^. The in-phase quantum ffuctuations have 
been over-compensated in this case: 



At 



(37) 



This indicates, that the reversal of the diffusion effect 
has been achieved by answering the fluctuations of the 
in-phase quadrature with an opposite field of double 
strength, effectively reversing the sign of that field com- 
ponent, while the unknown fluctuations in the out-of- 
phase quadrature are unmodified. 

Another interesting scenario is the stabilization of Sz = 
0, because it corresponds to the simple minded compen- 
sation of quantum fluctuations by chosing a feedback field 
with the negative amplitude of the quadrature measured 
in the homodyne detection. The effective diffusion step 
now reads 




Tt- 



An 



eff. 



X 

— S^Sy 



(38) 



feedback 



As shown in figure |^, this diffusion law has rotational 
symmetry around the Sx axis. The ratio of Sy/sz is a 
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constant. The diffusion is always directed towards one 
of the two stabiUzed poles with Sx = ±1- The absolute 
value of the diffusion constant for the stabilization of Sz = 
is thus given by 



r(i 



(39) 



In the general case of a stabilization of — cos 6, 
where 9 is the angle between the stabilized Bloch vector 
and the Sz axis, the diffusion step is 



6sx 



= Vrv— s 

a. 



eff. 




(40) 



and the diffusion constant is given by 



r(i 



2cos6is;, +cos^6'(4 + s^)). (41) 



By varying the feedback it is therefore possible to sup- 
press the diffusive dynamics for an arbitrary state in the 
Sx,Sz plane. 



C. Effects of the drift terms on timescales of l/F 

Although the feedback described above can completely 
suppress the random walk dynamics induced by the quan- 
tum fluctuations on a timescale of r for an arbitrary sys- 
tem state, it is necessary to consider the effects of the 
drift terms if stability on longer time scales is to be ob- 
tained as well. The complete diffusion step, including a 
feedback field of /(An) = — (1-1- cos 6*) An/2 | a |, is given 
by 



6Sy 

Ss, 








(42) 



where Ssq denotes the diffusion step without feedback, as 
given by equation (|20|). The drift is consequently modi- 
fied by the non-zero average of An in the feedback term. 
The resulting drift is the sum of the drift without feed- 
back and the feedback drift term. It may be written as 
the sum of a rotation around the Sx axis and a rotation 
around the s„ axis, 



(SSx) 

{SSy) 

(SSz) 



^(2 + S.)Sy 



Tt 



^ 




Sz 


) 








^ Sz 


9)sx 










(43) 



For the stabilized state with = sin 0, Sy = 0,Sz = cos 9, 
the effective drift is equivalent to a Rabi rotation around 
the Sy axis of (rr/2) sin0cos6' per time interval r. This 
effect can be compensated by a constant driving field. 
The field strength /o necessary for this purpose is given 
by 



/o 



Tt 



(44) 



This driving field ensures that the stabilized state is a 
stationary solution of the drift-diffusion dynamics gen- 
erated by the feedback setup. However, the drift terms 
may also cause problems if they amplify small deviations 
from this stabilized state. 

The stability analysis near Sx = siii9,Sy = 0,Sz = 
cos 9 can be performed by linearizing the drift dynamics 
of small deviations. The deviation in the Sx,Sz plane is 
appropriately described by the angular variable e, such 
that Sz = cos{9 + e) and Sx — sm{9 + e). The linear 
stability analysis then shows that for small e, the drift 
< Se > per time interval r is 



(Se) 



Ft 



(45) 



Therefore the drift terms always stabilize the state with 
suppressed fluctuations against rotations in the Sx , Sz 
plane. In the Sy,Sz plane, the situation is different how- 
ever. The linearized drift dynamics of small Sy is de- 
scribed by 

Tt 

(Ssy) ^ —{2 + cos9)cos9sy. (46) 

This indicates stability for all cos 9 < 0, i.e. any state 
which has a negative Sz component. If states of higher 
excitation are to be stabilized, any small deviation from 
Sy = in the initial preparation is amplified exponen- 
tially on a timescale of l/F. This indicates that long 
term stability of such inverted states cannot be achieved. 

An interesting aspect of the quantum trajectories as- 
sociated with the measurement protocols obtained from 
the homodyne detection is revealed in this critical prob- 
lem of quantum control. If the initial state is only known 
within an error margin, the trajectories may diverge and 
amplify this error margin even though the trajectories 
are fully deterministic functions of the initial state. Ap- 
parently, there is an aspect of deterministic chaos in this 
small quantum system when its dynamics is analyzed us- 
ing homodyne detection. It seems that the fully polarized 
states, Sx — ±1, Sy — Sz = 0, are ideal for stabilization, 
since they can be observed in the measurement results. If 
the system drifts away from its stabilized state, this can 
be observed in the homodyne detection and may be cor- 
rected, either by applying static fields to shift the phase 
of the atomic system or by shifting the phase of the local 
oscillator. In the following we will investigate the exam- 
ple of feedback stabilization for these maximally coherent 
states. 
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D. Stabilization of the dipole eigenstates 

For cos 6* — 0, the stabilized states are Sx = ±1, the 
eigenstates of the dipole component oscillating in phase 
with the local oscillator. Also, this case is special because 
it corresponds to the classical idea of feedback compensa- 
tion: An+Snext — 0. At the same time, the average of An 
is a measure of , indicating both the sign of the Bloch 
vector component stabilized and the success or failure of 
the stabilization attempt. The modified total diffusion 
step including feedback is given by 




(47) 



The compensation term has been split into one part 
which compensates the drift of the total diffusion step 
and the remaining term which modifies the effective drift. 
The nature of the diffusion term has already been dis- 
cussed in section III. The diffusion steps are directed 



towards the poles of the Bloch sphere with a diffusion 
constant of r(l — s^). The average drift is given by 




SxS^ 
SySz{2 + s^) 

-4{2 + s,)-sl 



(48) 



Figure H shows the drift vectors on the Bloch sphere. The 
change in always has the same sign as Sx, indicating 
that this drift will cause an increase in the dipole ex- 
pectation value Sx until an eigenvector with Sx = ±1 is 
reached. This implies that the feedback setup actually 
creates the stabilized coherent state regardless of the ini- 
tial state. 

If the initial state is the ground state, it is sufficient 
to examine the drift term for Sj, = 0, since neither the 
diffusion nor the drift will create an Sy ^ 0. The drift 
term is then given by 




Ft 







(49) 



This drift equation has four stationary solutions, two un- 
stable ones at Sz = ±1, the ground state and the excited 
state, and two stable ones at Sx = ±1, the dipole eigen- 
states. Consequently the system state will be drawn to- 
wards the closest one of the two dipole eigenstates as 
soon as the diffusion has moved the state away from the 
destabilized ground state. If the phase of the local os- 
cillator is changed, whether intentionally or accidentally, 
this does create an Sj, 7^ component. Since there will 
always be limits to the phase stability of the local oscil- 
lator, it is essential that the effects of such deviations are 
understood as well. For Sz — 0, the system reacts with a 
diffusion in the Sx, Sy plane and a drift reducing Sz- For 
small values of Sy, the system will relax exponentially to 
Sz = — 2Sy. This induces a drift in Sy corresponding to 



(50) 



This stabilization of = is very weak compared to 
the exponential stability in Sz = 0. The dynamics of the 
relaxation for small deviations is approximately given by 



At) 



1 



^4Tt + Sy{0y 



(51) 



Figure ^ shows a comparison of this relaxation with the 
exponential relaxation of Sz- This comparison clearly 
demonstrates that the relaxation of Sy may indeed be 
separated from the fast relaxation of Sz ■ Note that since 
the relaxation of Sy is very slow, its effect can not be 
separated from the diffusion. The diffusion constant of 
TSy is actually quite similar to the drift term of Sy. Con- 
sequently, the trajectory of the system dynamics will be 
more complicated and unpredictable than suggested by 
equation (^ij) and figure ^. 

On time scales longer than l/F, the expectation value 
of the dipole variable Sx can be observed in the averages 
of An. This may serve as a control of the stabilization 
process and possibly as an additional tool to produce the 
correct phase relation between the local oscillator and 
the atomic system. Note that it is fairly safe to inter- 
pret small deviations from = ±1 as phase mismatches 
between the local oscillator and the atomic system since 
the phase stability given by Sy is so much weaker than 
the excitation stability given by Sz- 



E. Implications for coupled and multi-level atomic 

systems 

The major part of this paper is concerned with the pos- 
sibility of observing and manipulating the dynamics of a 
simple atomic two state system interacting with the light 
field continuum. The theory of photon emission used may 
also be extended to coupled and multi-level systems, as 
shown in [ pO[ for a simple quantum beat scenario. One 
important aspect of the effects of homodyne detection of 
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light field emissions from larger electronic systems, such 
as networks of coupled quantum dots, is that the con- 
tribution of the system dynamics must be considered in 
more detail. In fact, the presence of different frequen- 
cies in the system dynamics may effectively remove the 
difference between homodyne and heterodyne detection. 
In quantum beat scenarios such as the ones considered 
in |l0| , [l^ , the homodyne detection could be performed 
on only one of the emission channels, with consequences 
for the probability of detecting an emission in the other 
channel. If the coupling strength of the transitions is 
very different [ p^ , the quantum fluctuations coupling to 
the fast transitions will induce phase fluctuations in the 
dynamics of the system. Such phase fluctuations in the 
quantum beats of a system could be measured by homo- 
dyne detection. The resulting correlation between the 
measurement protocols of the homodyne detection and 
the probability of an emission from a local part of the 
system can be predicted using quantum trajectory for- 
malisms such as the one presented in this paper. The 
effective projective measurement base of an eight port 
homodyne detection which was discussed in Q may also 
be applied in place of the simple balanced homodyne de- 
tection scheme used here. This would restore the symme- 
try between Sx and Sj,, increasing the number of possible 
stabilization scenarios. 



V. INTERPRETATION AND CONCLUSIONS 
A. Interpretation of the light-atom interaction 

In the usual photon detection measurement the quan- 
tum fluctuations of the vacuum field state seem to have 
no effect. This impression is a result of the particle pic- 
ture interpretation associated with the photon measure- 
ment. If homodyne detection is used instead, the infor- 
mation obtained is mainly information about the quan- 
tum fluctuations acting on the system, with only small 
contributions from the in-phase dipole component of the 
atomic system. While the information about the time of 
photon emission is lost completely, the evolution of the 
atomic dipole oscillations may be reconstructed from the 
measurement protocols. 

However, the Bloch vector does not react to the quan- 
tum fluctuations in the same way as it would respond to 
a classical driving field. The fact that no energy may be 
absorbed from the quantum fiuctuations of the light field 
vacuum requires that the ground state is not influenced 
by the fluctuations. At the same time, the correspond- 
ing excited state is much more sensitive to the fluctua- 
tions. This asymmetry required by energy conservation 
may also be understood in terms of a correlation of the 
atomic system and the field. The atomic ground state 
is really a dressed state in which the field and dipole 
fluctuations are correlated so as to preserve the state of 
lowest energy. Consequently, there is a similar correla- 



tion in the excited state, which enhances the interaction 
with the quantum fluctuations. Of course, this correla- 
tion extends only to the quantum noise. Therefore, the 
coherent fields used to control the atom interact in a dif- 
ferent manner and can never compensate the dissipative 
effects. This can only be achieved by manipulating the 
electromagnetic vacuum itself, for example by putting 
the atom into a microcavity. 

The peculiar nature of quantum states and quantum 
mechanical uncertainty is also apparent in the asymme- 
try of the diffusion step. If the quantum state is consid- 
ered to be an objectively real description of the atomic 
system it is difficult to explain the dependence of the dif- 
fusion step on the phase of the local oscillator used in 
the homodyne detection, especially since the measure- 
ment could be performed a long time after the emission 
and at an arbitrary distance. Thus this simple balanced 
homodyne detection scenario also highlights the cpiste- 
mological nature of the wavefunction and the resulting 
non-locality. 

The analysis of the contributions to the spontaneous 
emission process observed on time scales of 1 /F reveals a 
surprisingly complex structure of the process, with con- 
tributions from quasi classical dipole radiation and non 
linear effects of the fluctuations. The interpretation of 
these terms is far from complete and reveals new chal- 
lenges to our physical understanding presented by this 
fundamental quantum mechanical process. 

B. Conclusions 

We have presented a completely quantum mechanical 
theory of homodyne detection, simulating the time re- 
solved observation of electromagnetic emissions from an 
atomic system at a quantum efficiency of 100%. The 
theory includes a model of the unitary evolution of the 
system-field correlation and applies projective measure- 
ments to the results of this temporal evolution. The 
information obtained about the dissipative dynamics of 
the atomic system has been applied in feedback scenarios 
which demonstrate that the dissipative dynamics can be 
modified by feedback to create and stabilize excited or 
coherent states of the system. The dissipative nature of 
the measurement can not be compensated however, and 
excited states show an instability with regard to small 
errors in the initial preparation on time scales of l/F, 
the natural lifetime of the excited state. In conclusion, 
we have shown that homodyne detection can be a useful 
tool in attempts to observe and control individual quan- 
tum systems. 
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FIG. 5. Visualization of the effective diffusion step includ- 
ing a feedback stabilizing Sz = 0. The representation is as in 
figure 1. 
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FIG. 1. Visualization of the diffusion step on the Bloch 
sphere. The diffusion is represented by lines oriented parallel 
to the direction of the diffusion with a length proportional 
to the standard deviation of the diffusion step, a) shows the 
projection into the Sy,Sz plane and b) the projection into the 
Sx, Sz plane. 

FIG. 2. Total drift on the Bloch sphere projected into the 
Sy,Sz plane. Since the drift is rotationally symmetric around 
the Sz axis, no extra figure is given for the Sx,Sz plane. 

FIG. 3. Relative contributions to the exponential decay of 
the atomic excitation for j (s) ]= 1. 

FIG. 4. Visualization of the effective diffusion step includ- 
ing a feedback stabilizing Sz = +1. The representation is as 
in figure 1. 



FIG. 6. Effective drift on the Bloch sphere including a feed- 
back stabilizing Sz = 0. a) shows the projection into the Sy, Sz 
plane and b) the projection into the Sx,Sz plane. 

FIG. 7. Comparison of the approximate relaxation dynam- 
ics of Sz(0) = 0.1 and Sy{0) = 0.1 in the presence of a feedback 

stabilizing Sz = 0. Sy does not change much while the expo- 
nential relEixation effectively reduces Sz to zero. 
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